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To describe the compressible, turbulent flow in a labyrinth seal, a three volume 
bulk flow model is presented. The conservation equations for mass, momentum 
and energy are established in every control volume. A perturbation analysis is 
performed, yielding zeroth order equations for centric rotor position and first 
order equations describing the flow field for small rotor motions around the 
seal center. The equations are integrated numerically. From perturbation pressure, 
the forces on the shaft and the dynamic coefficients are calculated. 
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Forces on the shaft in x,y direction 

shaft displacements 

direct and cross-coupled stiffness 

direct and cross-coupled damping 

pressure 

axial, radial and tangential velocity 

temperature 

density 

axial, radial and tangential coordinate 
time 

length of labyrinth chamber 
depth of labyrinth chamber 
shaft radius 
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nominal clearance 
seal fin thickness 
total seal length 

eccentricity ratio, perturbation parameter 
radius of circular shaft orbit 
seal clearance 

rotational frequency of the shaft 
precession frequency of the shaft 
pressure before labyrinth 
pressure behind labyrinth 
swirl 

perfect gas constant 

specific heat ratio 

shear stress 

friction factor 

Reynolds number 

hydraulic diameter 

Mach number 

relative surface roughness 

entrance loss factor 

loss factor at seal fin entrance 

laminar viscosity 

imaginary unit 


SUBSCRIPTS 
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I, II, III 


zeroth order 
first order 
stator 
rotor 
axial 

circumferential 

fluid- 

chamber, upper wall 
chamber, side wall 
Sine 
Cosine 

seal entrance 
relative 

Indices for the different control volumes 
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INTRODUCTION 


Labyrinth seals, commonly used in turbomachinery to reduce the leakage loss, 
can have a strong influence on the dynamic behavior of a turbine or a compressor, 
especially for high performance machines. In some cases, the fluid forces gene- 
rated in labyrinths caused severe rotor instabilities. Therefore, the design engi- 
neer needs information concerning the possible influences that a labyrinth can 
have on the dynamics of a rotating system. 

For mathematical modeling, the seal forces can be described in terms of stiffness 
and damping parameters, the so-called rotordynamic seal coefficients (see Eq.l). 
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These coefficients are input parameters for rotordynamic investigations and have 
to be determined by fluid mechanics methods. In contrast to WEISER and 
NORDMANN (1987,1988,1989), where we used Finite- Differencing methods (FDM) 
to solve the conservation equations in conjunction with a turbulence model, 
Florjancic developed a more conventional three volume bulk flow method which 
takes advantage from FDM calculations and reduces the computational effort 
dramatically. In 1987, KLAUK has already developed a three volume model for 
the compressible flow in labyrinths, which is extended in WEISER (1989). In 
this approach, the energy equation and the equation of state for a perfect gas 
is used in addition to the previously existing model. 


MATHEMATICAL MODELING 

Many attempts have been made for the calculation of labyrinth coefficients, 
e.g. see the publications of KOSTYUK(1975) , CHILDS(1984) or WYSSMANNU984). 
Flow visualization experiments show, that for look-through labyrinths, the flow 
field can be divided in two characteristic regions: a vortex flow in the seal 
chamber and a jet flow region beneath the seal strip and the groove. 


Therefore, the calculation domain is divided into three control volumes (CV), 
see Fig.l: one in the seal chamber accounting for the vortex flow and two 
volumes representing the jet flow beneath chamber and seal fin. Fig.l also gives 
an impression of the bulk velocities in the different volumes. 
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Control Volumes Oth Order Bulk Velocities 



Fig. 1: 3-volume bulk flow model 

For mathematical description, we use the momentum, continuity and energy 
equation in every control volume. Paying attention to the following assumptions: 

- compressible medium (p / const.) 

- constant laminar viscosity p 

- modelling of turbulence by wall shear stresses 

- no shear stresses within the fluid (exception: jet friction shear stress 
in the contact region of CV II and CV III) 

- no radial shear stresses 

- first order radial velocities only in CV II and III 
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we obtain the following set of equations to describe the turbulent, compressible 
seal flow (equations shown here for teeth-on-stator seal): 
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As already mentioned, the modeling of flow turbulence is performed using semi- 
empirical wall shear stress formulations which is the usual way in bulk flow 
analysis, see for example CHILDS(1984), SCHARRER(1987). The fluid shear stress 
in the contact region of CV II and CV III is described in the same manner as 
WYSSMANN0984) suggested. All shear stress components can be written in 
a generalized form: 


*lj = ^ P *1 V j V rel 


with 
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im 
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and Re = — ^ , 

r/p 

V re i is the relative velocity, V. the component of the resulting velocity vector 
which is in the same direction as the shear stress component to be calculated. 
Index i indicates the appropiate shear stress (s: stator-, r: rotor etc.) and index 
j gives the direction (z: axial, <p: tangential). With these definitions, we obtain 
the following relations: 
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The parameters introduced in the shear stress formulations according to MOODY 
(1944) are 

n=0. 001375 
m=0. 33333 

C l =1 ° 4 

c 2 =io 6 

For the fluid shear stress, m, and C 2 are zero; n is set to 3.762 (according 
to WYSSMANN). To account for the real flow situation in the labyrinth, n is 
multiplied with the factor P , which depends on geometry and performance data. 


BOUNDARY CONDITIONS 

For the solution of the differential equations, boundary conditions have to be 
specified. To calculate the pressure, density and axial velocity at the seal en- 
trance, we use the same relations as NELSON(1985) , calculating the entrance 
values depending on the Mach number: 
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At the entrance of every seal fin, there occurs a pressure drop due to the 
acceleration of the fluid. Therefore, the equations 7 are used in a modified form 
with a constant loss factor 
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The model parameters 3, 3 v and are determined with the help of the Finite- 
Difference programs of NORDMANN and WEISER0987, 1988, 1989). 


PERTURBATION ANALYSIS 

In order to determine the seal forces, we assume that the rotor performs only 
small motions around the seal center. Then we can describe the dependent 
variables with a perturbations series, truncated after the linear term: 


u = u Q + e Uj 


P = Pq + epj 


h = C r + e hj 


v = v Q + evj 


T= T 0 + e Tj 
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w = w Q + e w ( 


P = P 0 + e Pi 


As perturbation parameter, we choose the relative eccentricity e, which is the 
ratio of the shaft eccentricity and the nominal clearance: 


e 



( 10 ) 


Inserting these expressions into the governing equations, separating terms without 
e and with e 1 yields two sets of equations: the zeroth order relations describe 
the flow field for a centric rotor position and the first order terms stand for 
the flow situation at small eccentric shaft motion. 
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ELIMINATION OF TEMPORAL AND CIRCUMFERENTIAL DERIVATIVES 


While the zeroth order equations depend only on the axial coordinate z, the 
first order set contains also derivatives with respect to time t, tangential and 
radial coordinates. Assuming that the shaft moves on a small circular orbit 
(see Fig. 2), the seal clearance can be written as 

h = h Q + ehj = C r - X(t) cos<p - Y(t) sincp (11) 



Fig. 2: Circular shaft orbit 


For the dependent variables, we prescribe a corresponding solution: 
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Now, the circumferential derivatives can be calculated analytically. Separating 
the equations with respect to sine and cosine terms and introducing complex 
variables: 
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we obtain a now complex system of first order equations. 

As already mentioned, we assume a circular shaft orbit, which gives for the 
time dependence of the clearance 

^r 0 e int (14) 


Again, similar expressions for the variables are established: 
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Thereby, the derivatives with respect to time can be expressed analytically. Fi- 
nally, the radial derivatives in CV II and III are modeled with 
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After having performed these steps, the first order equations do only depend 
on z. The solution starts with a numerical integration of zeroth order equations 
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with respect to the boundary conditions. The first order solution using the per- 
turbed boundary conditions yields the pressure distribution in the seal for small 
rotor motions. The forces are obtained by a pressure integration: 
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Because we have to determine 4 parameters (K, k, D, d), the calculation of first 
order equations has to be done twice for two different rotor precession fre- 
quencies (Q=0 and 0=w). Then the rotordynamic coefficients can be calculated: 
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COMPARISON TO EXPERIMENTS AND OTHER BULK-FLOW-THEORIES 


1. Example 

BENCKERT(1980) investigated experimentally a three chamber look-through laby- 


rinth. Due to the design of his test rig 
data of the labyrinth are 

R : 150 mm 
L : 7.75 mm 
B : 5.5 mm 
S : 0.25 mm 
C r : 0.5 mm 
p •. 1.425 bar 

~a 

p b : 0.953bar 


he could only measure stiffnesses. The 
T : 300 K 

a 

working fluid: air 
number of chambers : 3 

No shaft rotation ! 
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Fig. 3 and 4 show the comparison to the experimental data for different preswirl 
conditions. While the 3-volume approach predicts a positive direct stiffness 
according to the measurements, the theory of CHILDS0984) gives the opposite 
sign for K. The agreement between experiment and 3-volume model is good, 
especially showing that the model presented is capable to ensure sufficient 
accuracy in the calculation of dynamic seal parameters even for short labyrinth 
seals where the approaches of CHILDS or WYSSMANN have some problems. 


Comparison of measurement, bulk theory of CHILDS 
and 3-volume model 



— Experiments 
“ BENCKERT 

" Theory of 
CHILDS 

x 3— volume 
model 


preswirl [m/s] 


Fig. 3: Direct stiffness K 

Comparison of measurement, bulk theory of CHILDS 
and 3- volume model 



Experiments 

D BENCKERT 

TiT" Theory of 
CHILDS 

T-:~3- Volume 
model 


preswirl [m/s] 


Fig. 4: Cross-coupled stiffness k 
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2. Example 


In a paper published 1986, WYSSMANN shows a comparison of measured data 
for a 15-chamber look-through labyrinth. The experiments were carried out by 
CHILDS(1984). WYSSMANN compares the measurements to his 2-volume approach. 
For a 16-strip seal with fins on stator, the results of the experiments, of WYSS- 
MANN and the 3-volume model are shown in Fig. 5-7. 


Comparison of measurement, bulk theory of WYSSMANN 
and 3-volume model 



Experiments 

— CHILDS 
□ 

Theory of 

+ WYSSMANN 

3-volume 

model 


preswirl Lm/sJ 


Fig. 5: Direct stiffness 


Comparison of measurements, bulk theory of WYSSMANN 
and 3-volume model 



preswirl fm/sj 


Fig. 6: Cross-coupled stiffness 
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Comparison of measurement, bulk theory of WYSSMANN 
and 3-volume model 
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3-volume 

model 


preswirl [m/s] 


Fig. 7: Direct damping 

In contrast to the 2-volume theory, the model presented here predicts the direct 
stiffness more accurate. 


CONCLUSION 

A new bulk flow model has been presented for the calculation of the rotor- 
dynamic seal parameters. The theory uses a three volume approach, and conser- 
vation equations for mass, momentum and energy are set up for every control 
volume. A perturbation analysis yields zeroth and first order equations. The 
integration of first order pressure perturbation results in the desired dynamic 
coefficients. 

The comparison to experiments shows good agreement; compared to other bulk 
flow models, the present theory allows a much better calculation of the direct 
stiffness. This is mostly due to the determination of the model parameters with 
the help of the FDM programs, as a comparison to an earlier 3-volume model 
of KLAUK, which was developed in 1987 and did not have the benefit of using 
FDM results, shows. 
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